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g : ABSTRACT 

p ,1 We discuss the neutrino signature of a relativistic jet propagating through 

' a stellar envelope, a scenario realized in the collapsar model for Gamma Ray 

'sj" ' Bursts (GRBs). It is shown that the dramatic slowing of the jet deep within the 

^ ' star is accompanied by inelastic neutron-nucleon collisions and the conversion of 

^ a substantial fraction of the jet kinetic energy to neutrinos. These neutrinos have 

(X) ' observed energies in the range two to tens of GeV and an estimated detection 

^ ■ rate comparable to or larger than the detection rate of GeV neutrinos from other 

Tj- ■ GRB-related processes. The time delay between the arrival of these neutrinos 

^ ' and the GRB photons is tens of seconds. An observation of this delay would 

O ■ provide an indication that the GRB jet originated in a massive star. 

Subject headings: gamma rays: bursts — neutrinos 
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1. INTRODUCTION 

In this paper we discuss neutrino production in collapsars (Woosley 1993; MacFadyen 
& Woosley 1999). Collapsars, or failed supernovae, are one of two or a few currently favored 
explanations for the origin of long duration GRBs. In collapsars, the relativistic jet which 
shocks and ultimately gives rise to observed GRB photons must first pass through a massive 
stellar mantle. Here it is shown that this passage is associated with an interesting neutrino 
signal. 

A number of GRB-related neutrino sources have been previously investigated. These 
sources can be roughly divided according to whether they involve pion production via pho- 
tomeson interactions or via strong interactions. Neutrinos resulting from photomeson in- 
teractions typically have very high energies (> lO^TeV) and are possibly the most easily 
detectable in planned km-scale neutrino detectors such as ICECUBE (Halzen 1999). Wax- 
man & Bahcall (1997) discussed photomeson production of neutrinos in collisionless GRB 
shocks and argued that detection of these neutrinos would give information on the spectra of 
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photons and protons within the shocks, as well as place stringent constraints on fundamental 
neutrino properties. 

Neutrinos from inelastic nuclear coUisions typically have energies of the order of 10 GeV. 
Unless the phototube density in ICECUBE and other next generation neutrino telescopes is 
higher than currently planned, these detectors will have a small effective detector area for 
such low energy neutrinos and will not be sensitive to them (Karle 2002). However, there 
is a growing indication that, like neutrinos from p — 7 interactions, these neutrinos may be 
an important clue to conditions in and around the ORB central engine. Bahcall & Meszaros 
(2000) discussed the neutrino signal from inelastic collisions occurring during the dynamic 
decoupling of neutrons in the acceleration stage of the evolution of a GRB fireball. These 
neutrinos offer an indication of the composition and Lorentz factor of the GRB jet. Neutron 
diffusion in GRB environments was considered in a more general context by Meszaros & Rees 
(2000). Those authors showed that internal shocks at radii r ~ 3 • lO^^cm and transverse 
diffusion of neutrons into the jet at much larger radii can lead to an appreciable neutrino 
signal. The possibility of using neutrinos to distinguish between collapsars and supranovae 
(Vietri & Stella 1998, 1999) was discussed by Guetta & Granot (2002). The scenario 
envisaged by Guetta and Granot is unique in that it predicts TeV neutrinos from nuclear 
processes. This is because of the presence in their model of protons accelerated to very high 
energies in a pulsar wind bubble. 

The multi-GeV collapsar neutrino signal we study here is related to the mechanism 
giving rise to the variation in Lorentz factor of the outgoing jet. As with any GRB model 
invoking internal shocks, the jet variability must be substantial in order to explain observed 
temporal variability in GRBs, as well as to allow efficient conversion of kinetic energy to 
electromagnetic energy (see Piran (1999) for a general overview). In most GRB models, 
the jet variability is related to a variability in the mass ablation rate or luminosity of the 
central source driving the outflow. By contrast, the jet variability in the collapsar can arise 
from instabilities along the jet-stellar wall interface (Zhang, Woosley, & MacFadyen 2003). 
These instabilities operate at radii large compared to the size of the central black hole. 
Consequently, the instabilities act to slow a jet which has already had ample time to reach 
ultra-relativistic velocities. We point out that when free neutrons are present in the jet, this 
slowing is not elastic and is accompanied by a sizeable neutrino flux. 

In the next section we discuss details of how the slowing of the neutrons in the jet leads 
to the conversion of jet energy to potentially observable neutrinos. Estimates of neutrino 
detection rates and energies are presented. We will show that the neutrino signal of the 
slowing of the jet can originate from radii some two orders of magnitude smaller than the 
signal of processes discussed by Meszaros & Rees (2000). Also, the neutrino production we 
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discuss here is relatively insensitive to the bulk Lorentz factor of the outgoing jet, unlike 
other proposed sources of multi-GeV GRB neutrinos. 

2. Pion and Neutrino Production Deep Within the Collapsar 

The mechanism giving rise to the variability in the outgoing jet can be crudely repre- 
sented as the closing of a dense door composed of stellar material at some point along the jet 
axis. This results in a forward shock propagating into and accelerating the stellar material, 
and a reverse shock propagating backwards and slowing the jet. Because the stellar material 
blocking the passage of the jet is typically very dense relative to the jet material, the for- 
ward shock will only be mildly relativistic and uninteresting as far as a detectable neutrino 
signature is concerned. The top panel of Fig. 1 illustrates the shocking and structure of the 
outgoing jet. 

As particles in the jet traverse the reverse shock, they are decelerated from a Lorentz 
factor Tj^^ to a Lorentz factor F^j^'' < F^'^''. Here the superscript (L) denotes a quantity as 
measured in the rest frame of the collapsing star. The relative Lorentz factor across the 
shock is Trei ~ r^-^'*/2rit'*- Numerical simulations of relativistic jets propagating in stars 
indicate that the typical relative Lorentz factor is substantial. Zhang et al. studied the 
propagation of a relativistic jet in a 15Mq star for 3 different initial jet conditions. Typical 
values of Tj-ei for their simulations are given in Table 1. The simulations of Zhang et al., as 
well as observations of many peaks in observed GRB lightcurves, indicate that a substantial 
fraction (half or more) of the jet undergoes a slowing via relativistic shocks. 

The gross properties of the reverse shock slowing the jet are given by the Rankine- 
Hugoniot jump conditions. In particular, the proper baryon number density in the shocked 
jet (i.e. the baryon number density as measured in a frame comoving with the shocked fluid) 
is Ugh ~ ATreinj, where Uj is the proper baryon number density in the unshocked jet. For radii 
larger than ^ 10^F^-^''cm, which is the region of interest in the present work, the outgoing jet 
is coasting rather than accelerating (Piran, Shemi, & Narayan 1993). In the coasting regime, 
Tij is determined by baryon number conservation as nj = L/Airr^mNC^iXf'^)'^- Observations 
of GRBs suggest that the isotropic equivalent luminosity L is of order L > lO^^erg/sec (Frail 
et al. 2001), and that typically vf^ > 200 (Lithwick & Sari 2001). The energy density in the 
shocked jet fluid is Ugh ~ Fre/ris/im^vc^, provided that the specific enthalpy in the unshocked 
jet is small (i.e. that the jet is in the coasting regime). In the middle panel of Fig. 1 we 
show the evolution of rij and F^^^ in the jet. 

Pion production and associated neutrino production in the reverse shock depend in detail 
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on the processes mediating the shock. It is convenient in discussing the shock structure to 
consider separately those particles which interact cicctromagnetically (protons-e^-photons) 
and the neutrons. Neutrons can only be slowed by strong collisions with other nucleons. 
Because the length scales associated with electromagnetic processes are small compared to 
the neutron mean free path, an upper limit to the width of the shock can be obtained by 
considering the case where strong nucleon-nucleon collisions alone mediate the shock. Such 
coUisional shocks have been studied by Mott-Smith (1951) in the non-relativistic case and 
by Chaphne & Weaver (1979) and Cercignani & Majorana (1988) in the relativistic case. 

Wc take as an upper limit to the shock width the collisional width 5c ~ l.'^/Vreinj'yNN 
found from the Monte-Carlo simulations of Chapline & Weaver (1979). Here (Jnn is the 
strong scattering cross section and the shock width is as measured in a frame comoving with 
the shock. For center of mass energies a few hundred MeV above the threshold for pion 
production, (Jnn ~ 4 • 10~^^cm^, and approximately 80% of this cross section goes towards 
inelastic production of pions (see Hagiwara et al. (2002) for a compilation of strong cross 
sections). For this work, the important properties of 5c are I) that because the shock width 
is at most a few collisional mean free paths, the neutron-nucleon collisions slowing the shock 
are hard and produce pions, and II) because 5c oc l/nshi^NN-, the number of inelastic collisions 
suffered by a typical neutron in traversing the shock is approximately independent of F^e/- 
A close-up view of inelastic neutron-nucleon collisions and associated pion production in the 
reverse shock is shown in the bottom panel of Fig. 1. 

Protons in the shock will be slowed by electromagnetic, rather than strong, processes. 
Deep within the star, where the shocks considered in this paper are occurring, the plasma 
is optically thick and the thermalization time for radiation is short. Consequently, we ap- 
proximate the protons as being in thermal equilibrium with the background e^/7 plasma^. 
The temperature of the protons in this case is small, kT <^ mpC^, and the protons are 
non-relativistic in the shocked plasma frame. 

Our assumption about the proton distribution function is in accord with results from 
studies of shock propagation in type II SNe (e.g. Ensman & Burrows (1992)). However, 
Meszaros & Waxman (2001) have argued that protons shocked in mildly relativistic collisions 
at large radii inside the coUapsar will be driven to a power law distribution, but that protons 
shocked in the strongly relativistic colhsion between the jet head and the stellar mantle may 
or may not have a power law distribution. The question of the precise radius (presumably 



"'^The opposite limiting case for the proton distribution is that where the protons are Fermi- accelerated 
and achieve a power law distribution. In this case, the maximum proton energy (typically much larger than 
the proton rest mass) is that for which the acceleration and cooling times are equal (Waxman 1995). 
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larger than lO^^cm) at which shocks are coUisionless rather than radiative is not addressed 
here. As we are most interested in processes at r < lO^-'^cm, this issue is not so important 
for the present study. 

In the inelastic collisions slowing the neutrons, the total energy going into pion produc- 
tion in the center of mass frame of the coUision is ^ (l/2)(£'cm — Sm^^c^) (e.g. Mannheim 
& Schlickeiser (1994)). Here Ecm is the energy of the collision in the center of mass frame, 
and the above equation is valid provided that Ecm — 'imj^c^ is larger than a few times the 
pion rest mass energy. Very near the threshold for pion production, E^°^ ^ i^cm — 2mj^c^. 
The average energy of individual pions created in the collision is E^°^ /^^ where ^ is the mul- 
tiplicity of produced pions. As measured in a frame at rest with respect to the star, the 
energy of pions produced in a neutron-nucleon collision is 

= T^mEr' = (^) (rf ) + r^) - 2r.^) . (i) 

Here we have made the assumption that the typical collision occurs between a nuclcon with 
Lorentz factor F^^'' and a nucleon with Lorentz factor T^j^'*. In Eq. 1, Tcm is the Lorentz 
factor of the center of mass frame as measured in the lab frame. In terms of the Lorentz 
factors of the unshocked and shocked portions of the jet, Vcm ~ ^'Tp^^T^ if Vgh is larger 

than two or three, and Vcm ~ \J r^^^ / 2 if Tgh ~ 1- The observed energy of an individual pion 
created in a collision will again be a factor of 1/.^ smaller than the estimate in Eq. 1. For 
example, for F^'^'' = 100, eI^^ pa 10, 12, and 7 GeV, for Tgh = 30, 10, and 1 respectively. 

The neutrino signal of the inelastic collisions depends on how the spectrum of pions 
produced in these collisions is modified by cooling processes. Pion cooling may proceed 
through inelastic strong scatterings, or via electromagnetic interactions - coulomb scattering 
off the background electrons, inverse compton scattering, and synchrotron emission. 

Cooling via strong scatterings is simply characterized by the vr— nuclcon interaction 
time. We can roughly approximate tt — N collisions as occurring in the shocked region (as 
opposed to the region of width 6c characterizing the shock). In this case the tt— nucleon 
interaction time is ^ 

Tn-N = , (2) 

CO-Nnnsh 

where a^v^ is the nuclcon-pion scattering cross section. For center of mass energies a few 
hundred MeV larger than the pion rest mass, (Jntt ~ 2.5 ■ lO^^^cm^ (Hagiwara et al. 2002). 
As with NN coUisions, the bulk of this cross section goes toward inelastic tt production. 
Though inelastic ttN collisions may roughly conserve the total energy in pions, they tend to 
reduce the average pion energy, making detection of the product neutrinos more difficult. 
Note that because (Tntt and (Jnn are comparable, a more accurate expression for r^-jv should 



-6- 



involve an effective nucleon number density witfiin tfie region of width Sc (rather than the Ush 
appearing in Eq. 2). Because the number density within the shock is smaller than Ugh, Eq- 
2 represents a lower bound on t^^-n (or an upper limit on the influence of ttN scatterings). 

Pion energy loss via coulomb scattering off the background electrons is given by the 
Bethe-Bloch formula for ionization losses (see Groom et al. (2000)), with the ionization 
potential for completely ionized material equal to fkUp. Here ujp = ^J ^.Tm^e^ j mf. is the 
electron plasma frequency, and is the electron number density. Neglecting an unimportant 
(for our applications) logarithmic dependence on rig, pion energy loss via coulomb scattering 
is given by dE/dx (2/3) ■ 10~^^neGeV/cm. The timescale for energy loss via tt — e~ 
scatterings is then Tt^-^- ~ (r^/2)(nsh/'n.e)T7r_Ar, with V^^ the pion Lorentz factor as measured 
in the shocked fluid frame. When strong scatterings are inefficient at cooling the pions, the 
temperature in the shocked plasma is generally too low for efficient pair production. In 
this case rie ^ Ugh and tt — e~ interactions are at most only as important as strong scatterings. 
Pion energy loss via bremsstrahlung and pair production is not important for pion energies 
less than about 100 GeV (Groom et al. 2000). 

Pions lose energy to the background photon and magnetic fields at a rate {dET^/ dt)^^^^ 
(cr^)(me/m,r)^crrC^em, where (Tt is the Thomson cross section, m,r is the pion mass, and Uem 
is the energy density in the photon and magnetic fields. A cooling time be defined 



over cooling via strong interactions when 



as T^ — 'y — E^{dE^/dt)-^^. Coohngof pions via these electromagnetic interactions dominates 



< 1- (3) 



400/ KushmNC^J V 400 



Here we have made use of the fact that Uem ~ Ush for relativistic shocks in an optically 
thick flow. Only in some of the more extreme cases that we consider here is the above 
inequality not satisfied. By contrast, for shocks at large radii (r ~ lO^^cm) when the protons 
are accelerated to high energies via plasma processes, typical pions produced in photomeson 
interactions have Lorentz factors of ~ 10^, and synchrotron cooling is more important than 
strong scattering. 

In the shocked plasma rest frame the lifetime of a charged pion is t^^-d = 2.6 ■ lO^^FTrS. 
The fraction of pions that decay before suffering a strong scattering is I/Ittt-d/ttt-n + !)■ 
The pion Lorentz factor relative to the shocked plasma for which r^-z? = t^^-n is 

_ (r/10^°cm)^(rf)/100)^ 

J- 17,0001 = I ■ i-O . 1,4 j 

Pions with a Lorentz factor greater than r^^^cooi will on average scatter before decaying, and 
pions with a Lorentz factor smaller than cooi will on average decay before scattering. Eq. 
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4 indicates that for r > lO^cm, pions produced in mildly relativistic collisions (Trei ~ 2) 
will decay before cooling. Conversely, for r ~ lO^cm, pions produced in strongly relativistic 
shocks {Trei ^ 10, as in case JA in the calculations of Zhang et al.) will cool before decaying. 
Neutrinos from the decay of cooled pions will typically have observed energies < IGeV, and 
will be difficult to detect. For r > lO^^cm, Eq. 4 indicates that pions produced even in very 
relativistic shocks will decay before cooling. 

In Fig. 2 we show the ratios of the timescales for the pion cooling processes discussed 
above to the timescale for pion decay. For this figure we have approximated the electron 
and nucleon number densities in the post-shock region as being equal and we have also 
made the assumption that Trei = F^. As can be seen from Fig. 2, for Trei = 2 and 
r > 1.5 • 109v'-^52(100/F^-^'')cm, pion cooling processes are unimportant (e.g. the timescale 
for all coohng processes is more than twice the decay timescale and less than 30% of the 
pion energy is lost). For Frei — 10, pion cooling processes are unimportant for r > 7 • 



On average, neutrinos carry away ~ 2/3 of the pion energy. The typical energies of 
the three neutrinos the pion (and daughter muon) decays to are in the range 30-50MeV 
as measured in the pion rest frame. This implies observed neutrino energies in the range 

~ [1/5 — l/3](£^i^''/^)(l + z)~^GeV, with z the redshift at which the burst occurs. In Table 
1, typical neutrino energies for the different simulations of Zhang et al. are shown. Note 
that the observed neutrino energy scales linearly with jet Lorentz factor. For example, for 
F^^^ = 400, the expected neutrino energies are a factor of approximately four higher than 
shown in Table 1. 

The fraction of jet kinetic energy lost to neutrinos is approximately 



caying. Here the electron fraction Ye is the ratio of protons to baryons in the fiow. Recent 
studies indicate that the outgoing jet is likely neutron rich, with plausible values of Y^. in 
the range 1/20 < i; < 1/2 (Pruet, Woosley, & Hoffman 2002; Beloborodov 2002). In Eq. 
5, one factor of 2/3 approximately accounts for the fraction of charged pions (as opposed to 
7r°'s) produced, while the other factor of 2/3 accounts for the fraction of the charged pion 
energy going to neutrinos. The estimate in Eq. 5 assumes that after one inelastic collision 
the decelerating neutron does not undergo subsequent inelastic collisions. For this reason, 
Eq. 5 represents a lower bound to the fraction of jet energy lost to neutrinos. The true 
fraction could be as much as a factor of ~ 2 higher, depending on the jet composition and 




cm. 




(5) 



for collisions ocurring at radii large enough that the charged pions do not cool before de- 
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details of the shocking. For the jet dynamics calculated by Zhang et al., estimates of the 
fraction of jet energy going to neutrinos lie in the range [0.1 — 0.2] (1 — Yg)- 

To facilitate comparison with studies of other GRB related neutrino production pro- 
cesses, we closely follow the analysis and notation of Bahcall & Meszaros (2000) and Meszaros 
& Rees (2000) in estimating the neutrino detection rate. In the notation of those authors, the 
average neutrino detection rate of events at rcdshift z is ^ (Nt/A7rD'')RbNna,,. Here Nt 
is the number of protons in the terrestrial detector, D is the proper distance out to rcdshift 
z, Rb is the GRB rate within a Hubble radius, and Nn is the isotropic equivalant number 
of neutrons undergoing inelastic scattering. The neutrino detection cross section averaged 
over neutrinos and anti-neutrinos is cTj^p ~ 0.5 • 10~^^{'^Ei,/lGeV){l + z)~^cm^ (Gaisser 
1990), where ~ (2/3)M^^ is the total energy of all neutrinos resulting from an in- 

elastic neutron-nucleon collision. Using Eq. 1 and adopting an Einstein-de Sitter cosmology 
with Hubble constant H — dbhQ^km/s/Mpc gives 

. .... - n, (S) m m (TTfe^b)^- 

(6) 

Here £^53 is the isotropic equivalant jet energy in units of lO^^erg, and f shock is the fraction 
of the jet undergoing a relativistic slowing at radii large enough the produced pions do not 
cool and decay to unobservable neutrinos (r > lO^cm for mildly relativistic collisions and 
r > lO^'^cm for ultra- relativistic collisions). For illustration, the models of Zhang et al. 
predict neutrino detection rates of ~ 21/year (JA), ~ 7/year (JB), and ~ 10/year (JC), 
for £"53 = 1, = 1/2, and fshock—^/'^- For model JB, where the relative Lorentz factors 
across the shock deep within the star are very similar to the relative Lorentz factors typical 
of internal shocks occurring at ~ lO^^cm, the expected detection rate is roughly the same as 
that calculated by Meszaros & Rees (2000) for neutrinos from internal shocks. For models 
JA and JC, in which the jet is more dramatically slowed, the expected detection rates are 
somewhat larger, approaching the sorts of detection rates expected for efficient sideways 
diffusion of neutrons into the jet. 

As can be seen from Eq. 6, R,y does not depend explicitly on the beaming angle of emit- 
ted neutrinos. However, as we discuss below, a temporal and directional correlation between 
detected neutrinos and the parent GRB is needed in order to distinguish GRB neutrinos 
from background atmospheric neutrinos. If the beaming angle of emitted neutrinos (D,^) is 
smaller than the beaming angle of GRB photons (Jl-y), then a parent GRB is always visible 
in coincidence with detected GRB neutrinos. If the opposite case holds {fli, > fl^y), then 
those neutrinos emitted outside of the photon beaming angle are effectively lost to observa- 
tion and the detection rate is effectively decreased by a factor of fl^/fl^. Two competing 
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processes affect the focusing of an initially wide jet as the jet propagates through the 

star (Aloy et al. 2000; Zhang, Woosley, & MacFadyen 2003), and the expansion of the high 
enthalpy jet as it breaks free of the star (Zhang, Woosley, & MacFadyen 2003). Jet focusing 
tends to increase Vtu/Vt^ because it makes the jet opening angle within the star larger than 
the jet opening angle at the stellar surface. Expansion of the high enthalpy jet tends to 
decrease Vly/Vl^. Present 2-d simulations of jet propagation do not definitively determine 
which process is dominant, though they are roughly consistent with ^ 1. 

It should be emphasized that Eq. 6 and the detection rates estimated for GeV neutrinos 
from GRBs neglect many of the experimental difficulties associated with detecting these 
neutrinos. As noted in the introduction, the detector phototube spacing must be small 
enough to allow detection of ~ lOGeV muons with ~ 50m pathlengths (e.g. see Halzen & 
Hooper (2002)). In addition, the angle of the incident neutrino relative to the GRB must 
be sufficiently well known to distinguish a GRB neutrino from the background atmospheric 
neutrinos. Rough estimates can be made by scaling detection rates for Supcr-Kamiokande 
to the larger detectors discussed here. In the search through the Supcr-Kamiokande data for 
GRB neutrinos, a high-energy neutrino {E,y > 200MeV) background rate of 2 ■ 10~^(20 sec)^^ 
was assumed (Fukuda et al. 2002). Estimating that roughly 1/5 of these have energies larger 
than a few GeV, approximating the volume of a km^ detector as being ~ 10^/5 times larger 
than the volume of Super-Kamiokande, and supposing the angle of the incident neutrino 
relative to the GRB position is known to within a cone of ~ 10° half-opening angle gives a 
background in a km^ detector of 0. 1(20 sec) With this background rate, a total sample 
of 1000 GRBs, and no other information, a detection of 20 neutrinos from 20 separate GRBs 
(out of a thousand) is not statistically significant, but the detection of two or three neutrinos 
from a single GRB is quite significant (see Table 1 in Fukuda et al. for detailed estimates). 
As a definite example of the potential of large neutrino detectors, the slowing of a relativistic 
jet in GRB 980425 {z = 0.0085 (Galama ct al. 1998)) should have given roughly 100 neutrino 
detections in a km^ detector. Even without directional information, such a detection would 
be significant. 



3. Summciry 

Numerical simulations of the propagation of a relativistic jet through a stellar envelope 
indicate that, deep within the star, instabilities along the jet-stellar wall interface slow the 
jet. This slowing gives rise to a variation in the Lorentz factor of an initially uniform jet. As 
we have shown, the slowing of neutrons within the jet also gives rise to an interesting and 
potentially detectable neutrino signal. 
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Neutrinos produced during the slowing of the jet can precede GRB photons by tens to 
several tens of seconds, depending on the jet dynamics and stellar progenitor. By contrast, 
neutrinos produced in a jet that is not subsequently dramatically slowed arrive within ~ 
0.3(r^/10-'^^cm)(100/r^-^^) seconds of the GRB photons, where is the radius where the jet 
shocks and produces observed photons. The observation of GeV neutrinos preceding a GRB 
flash by ~ 20 seconds would be an indication that a relativistic jet had been slowed while 
optically thick. In turn, this could imply that the jet had punched through a stellar envelope. 

Slowing of the jet within the star also has an impact on the jet composition. This 
is because neutron-neutron collisions, which dominate the slowing of neutrons if the jet is 
neutron rich, preferentially produce protons. Near threshold, for example, roughly 75% of 
the inelastic nn(pp) cross section goes toward the production of a p(n) (McGill et al. 1984). 
A jet born near the central black hole with a very low electron fraction will arrive at large 
radii with ~ 1/2. This will influence subsequent nucleosynthesis in the jet (Pruet, Guiles, 
& Fuller 2002; Lemoine 2002), as well as the production of neutrinos by inelastic nuclear 
processes at large radii. 

Apart from composition effects, though, the source of GeV neutrinos we have discussed 
arises independently of other proposed sources of GeV neutrinos in relativistic fireballs. For 
example, neutron diffusion in internal shocks at r > lO^^cm or the transverse diffusion of 
neutrons into the jet at larger radii (Mcszaros & Rees 2000), and inelastic nuclear collisions 
during neutron decoupling (Bahcall & Meszaros 2000), can still occur and produce neutrinos. 
In the most optimistic case where all of these various processes occur, the detection rate of 
neutrinos in a km-scale detector could be as high as 30 per year. 
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conditions in coUapsar jets and several useful suggestions from the referee John Beacom. 
This research has been supported by the DOE Program for Scientific Discovery through 
Advanced Gomputing (SciDAG; DE-FG02-01ER41176). This work was performed under the 
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Fig. 1. — Illustration of the jet shocking and shock structure. Panel a) shows the shocking 
of the jet propagating through the stellar envelope. In this paper we are concerned with 
shocks occuring between r ~ 10^ — lO^^cm. Panel b) shows the the baryon number density 
(n) and jet Lorentz factor (F) in the un-shocked (subscript j) and shocked (subscript sh) 
regions. Panel c) shows a close-up of the reverse shock (circled region in panel b). In panel 
c) the shock width is Sc- This is the typical length over which the neutron-nucleon collisions 
slowing the jet occur. Neutrons will eventually thermalize with the plasma in the shocked jet 
region. Depending on the details of the shocking, pions created in inelastic neutron-nucleon 
collisions may or may not thermalize before decaying. 
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Table 1. Estimates of properties of the neutrino signal from jets with different initial 

conditions propagating through a star. 



ModeP 


r t> 

rel 


E,{1 + zY 




JA 


~6-50 


2-3 


21 


JB 


~ 2 - 3 


5-8 


7 


JC 


^ 4 _ 9 


2-3 


10 



'^Models refer to the different models 
calculated in Zhang, Woosley, & Mac- 
Fadyen (2003) 

^Typical Lorentz factors across the 
shocks slowing the jet within the star 
from the simulations of Zhang et al. 

'^Our calculation of the typical range 
of energy in GeV of neutrinos created 
in the shock slowing the jet as seen by 
an observer at rest with respect to the 
star. 

'^Our estimate of the annual detec- 
tion rate of neutrinos with energies 
of a few GeV in a km-scale detector. 
In this estimate, z = 1, = 1/2, 
E = lO^^erg, and f shock = 1/2 are as- 
sumed. 
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Fig. 2. — Illustration of the ratio of the timescales governing the various coohng processes 
to the timescale for pion decay. Here a = (l/\/L^(r^-^^/100). The total cooling timescale 
is T;r_c = {t^-n + '''n-'y + '^w-e)"^- ^or the F^eZ = 2 curvcs, t-„_^/tt,_d is too large too appear 
on this graph, and the Ttt-jv/ttt-d curve exactly overlaps with the TT^^e/Tn-n curve. 



